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ABSTRACT 

Light travel time delays distort the apparent shapes of HII regions surrounding bright quasars during early 
stages of cosmic reionization. Individual HII regions may remain undetectable in forthcoming redshifted 21 cm 
experiments. However, the systematic deformation along the line of sight may be detectable statistically, either 
by stacking tomographic 21cm images of quasars identified, for example, by JWST, or as small-scale anisotropy 
in the three-dimensional 21cm power spectrum. Here we consider the detectability of this effect. The anisotropy 
is largest when HII regions are large and expand rapidly, and we find that if bright quasars contributed to the 
early stages of reionization, then they can produce significant anisotropy, on scales comparable to the typical sizes 
of HII regions of the bright quasars ( <S 30Mpc). The effect therefore cannot be ignored when analyzing future 
21cm power spectra on small scales. If 10 percent of the volume of the IGM at z — 10 is ionized by quasars with 
typical ionizing luminosity of S ^ 5 x 10 56 sec _1 , the distortions cause a <; 10 percent enhancement of the 21cm 
power spectrum in the radial (redshift) direction, relative to the transverse directions. The level of this anisotropy 
exceeds that due to redshift-space distortion, and has the opposite sign. We show that on-going experiments such 
as MWA should be able to detect this effect. A detection would reveal the presence of bright quasars, and shed 
light on the ionizing yield and age of the ionizing sources, and the distribution and small-scale clumping of neutral 
intergalactic gas in their vicinity. 

Subject headings: cosmology: theory - early Universe - galaxies: formation - high-redshift - evolution - 
quasars: spectrum 
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1. INTRODUCTION 

How and when the intergalactic medium (IGM) was reion- 
ized is one of the long outstanding questions in cosmology, 
likely holding many clues about the nature of the first gener- 
ation of light sources and the end of the cosmological "Dark 
Age" (see, e.g., Barkana & Loeb 2001 for a review). 

Observational break-throughs in recent years have revealed 
our first clues about this epoch. The Thomson scattering op- 
tical depth, r e = 0.09 ± 0.03 (Page et al. 2006; Spergel et al. 
2006), inferred from the polarization anisotropics in the cosmic 
microwave background (CMB) by the Wilkinson Microwave 
Anisotropy Probe (WMAP), suggests that reionization began at 
redshift z ~ 10. The detection of the so-called Gunn-Peterson 
(GP) troughs in the spectra of high-redshift quasars discovered 
in the Sloan Digital Sky Survey (White et al. 2003; Fan et al. 
2006), and more detailed analyses of these spectra (Mesinger 
& Haiman 2004; 2006), suggest that reionization is ending at 
z-6. 

Despite this progress, little is known about the nature of the 
first ionizing sources. The dearth of bright quasars at z ~ 6, and 
a tight limit on the contribution of lower-luminosity quasars 
to the unresolved soft X-ray background (Dijkstra, Haiman & 
Loeb 2004) requires that the ionizing background at z ~ 6 be 
dominated by stars, rather than quasars. On the other hand, 
luminous quasars may still contribute up to ~ 10% of the ion- 
izing background at z ~ 6 (e.g. Sbrinovsky & Wyithe 2007), 
and the contribution of the ionizing radiation from accreting 
quasar black holes may have been more significant at higher z, 
during the earliest stages of reionization (Madau et al. 2004; 
Ricotti & Ostriker 2004; Oh 2001; Venkatesan, Giroux & Shull 
2001). Overall, the present data allow a wide range of possible 



reionization histories, driven by different sources and various 
physical feedback mechanisms (e.g. Haiman & Holder 2003). 

One of the promising future probes that could constrain, es- 
pecially the early stages, of reionization, is based on the de- 
tection of the redshifted 21 -cm line of neutral hydrogen (see a 
recent review by Furlanetto, Oh & Briggs 2006). Forthcom- 
ing 21cm surveys are expected to deliver tomographic maps of 
the brightness temperature of the 21cm line, and detect fluctu- 
ations corresponding to the topology of large discrete HII re- 
gions. Numerous recent studies have addressed various aspects 
of the power spectrum of such brightness temperature fluctua- 
tions (see references in Furlanetto et al. 2006). 

In this paper, we examine one particular effect, arising from 
the distortion in the apparent shapes of HII regions due to the 
light-travel time delay between a light source and the ioniza- 
tion front it drives into the IGM. Since the distance from the 
Earth to various points on the surface of an HII region differ, 
different patches of an expanding HII region will be observed 
at different stages of their evolution, corresponding to different 
light travel times. As pointed out in Cen & Haiman (2000), 
the delay can be ignored along the line of sight to individual 
sources (for example, when computing Lyman a absorption 
spectrum of a point source; see also White et al. 2003). On 
the other hand, the delay modifies the apparent expansion of a 
resolved HII region in the transverse direction, as could be ob- 
served, for example, in 21cm studies (Wyithe & Loeb 2004c). 
The full equal-arrival time surface in three dimensions has been 
computed in simple models of individual spherical quasar HII 
regions by Yu (2005), and shown to be highly anisotropic for 
young and bright sources. 

In principle, the distortion could be measured directly in to- 
mographic images of individual HII regions (Wyithe, Loeb & 
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Barnes 2005), although the immediate next generation of 21cm 
instruments (such as PAST or MWA) are unlikely to be able to 
achieve the required S/N for such direct imaging. Since individ- 
ual HII regions are likely to be intrinsically highly anisotropic, 
such a measurement can only be performed statistically. In this 
paper, we examine and quantify the possibility of measuring 
the corresponding anisotropy in the three — dimensional power 
spectrum of the 21cm brightness temperature. This method 
does not require that individual HII regions be detectable at 
high S/N. 

Barkana & Loeb (2006) recently studied the anisotropy in 
the three-dimensional 21cm power spectrum that arises from 
finite light-travel time across larger, joint HII regions that have 
merged around galaxies. The growth of such super-bubbles is 
dictated by the collapse rate of cosmic structures, and it pro- 
duces an anisotropy which will be significant on larger scales 
(comparable to the sizes of the HII super-bubbles). In com- 
parison, the effect we describe here arises from the finite rate 
of growth of isolated HII regions around individual quasars. 
Finite-speed-of light effects during the lifetime of quasars can 
produce an appreciable additional smaller-scale anisotropy only 
in a likely narrow window of redshifts during reionization: on 
the one hand, bright quasars need to present, which occurs only 
at relatively late stages, but on the other hand, once the ion- 
ized bubbles around these bright quasars percolate significantly, 
the small-scale anisotropy will diminish. Conversely, a detec- 
tion of any small-scale anisotropy would reveal the presence 
of bright quasars, and shed light on the ionizing yield and age 
of the ionizing sources, and the distribution and small-scale 
clumping of neutral intergalactic gas in their vicinity. In this 
paper, we quantify the conditions under which a detection of 
this additional anisotropy, due to the presence of quasars, could 
be feasible with forthcoming experiments, with specifications 
similar to that proposed for MWA. 

The rest of this paper is organized as follows. In §[2] we dis- 
cuss the basic formalism to compute the 21cm brightness tem- 
perature fluctuations, including the anisotropy from the light 
travel time delay. In § [3] we present the expected level of 
anisotropy in a few simple models, and discuss how it is ex- 
pected to vary with various parameters. In §2] we discuss in- 
strumental noise and the corresponding limits on the detectabil- 
ity of the anisotropy. In § [5] we discuss our results and the im- 
plications of this work. Finally, in § [6] we summarize our main 
conclusions. Throughout this paper, we adopt the background 
cosmological parameters D, m = 0.29, £l\ = 0.71, = 0.047, 
Hq = 72 km s" 1 Mpc -1 , with a power spectrum normalization 
agffi = 0.8 and slope n s = 1, consistent with the values measured 
recently by the WMAP experiment (Spergel et al. 2007). Un- 
less noted otherwise, all lengths below are quoted in comoving 
units. 

2. BRIGHTNESS TEMPERATURE FLUCTUATIONS 

Forthcoming radio observations will measure the brightness 
temperature at the redshifted wavelength (1 +z)2 lcm. The bright- 
ness temperature is expected to be below that of the CMB be- 
fore the first luminous sources turn on, so that cosmic HI can 
only be seen in absorption. Once the first sources turn on and 
ionize and heat the IGM, the brightness temperature can rise 
above that of the CMB, so that the 21cm signal is seen in emis- 
sion (see Furlanetto et al. 2006 for a review and references). 
Here we are interested in the epoch of reionization, when the 
spin temperature T s ^> 7cmbr- In this limit, the brightness tem- 



perature T of the redshifted HI line is independent of T s , and 
is fully determined by the local density of neutral HI (see e.g. 
Field 1959; Madau, Meiksin & Rees 1997; Sethi 2005). In this 
section, we outline how we compute the anisotropic correlation 
function of the spatially varying brightness temperature. 

2.1. Basic Formalism for Temperature Fluctuations 

Observations will measure the brightness temperature T at a 
given frequency v and direction on the sky ft, but in a narrow 
frequency band, the signal is dominated by the emitting gas 
at the corresponding redshift z, or radial distance r (with the 
relation between r and observed frequency v obtained from the 
appropriate cosmological model). Given the telescope beam 
B(n,n'), which for radio interferometers will correspond to the 
synthesized beam, we obtain: 

T(n,r)= I dQ! [ dr%(/)Wr(r')ip(n'y)B(fi,n'), (1) 



where 7o(r) is the mean brightness temperature at the redshift 

Z = z(r), 

3 h\ 2 A 2l n h (z) 
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V=x H (r)[l + <5(r)], 



(3) 



with S and x H denoting the spatially fluctuating over-density and 
neutral hydrogen (HI) fraction (n m /[n m +n m J), respectively. Fi- 
nally, W r (r') is a window function, peaking at r, and represent- 
ing the bandwidth of the observations. 

Assuming for simplicity that 7o(r) does not evolve signifi- 
cantly within the bandwidth, it can be taken out of the integral. 
In this case, 

T(ft,r) = T (r) [ dfl' [ dr'W r (r')ip(n' ,r')B(n,n'), (4) 



where r corresponds (for example) to the center of the observed 
frequency band. The two-point correlation of the temperature 
fluctuation can then be written as 

(A7X"i,ri)A7Xn 2 ,r 2 )) = T ( ri )T Q (r 2 ) x 

xfdrj dr'J dQdSl'W n (r)Wr 2 (r')B(n,n 1 )B(n' ,n 2 ) 

[(iP(n u r)i;((n 2l ri))-(^) 2 ]. (5) 

Note that statistical homogeneity of the signal is assumed but 
not statistical isotropy, i.e. the correlation depends only on 
ri-r 2 , and the brackets denote averages over (say) ri. The ex- 
pression above also assumes that the neutral fraction does not 
evolve significantly over the length-scales of interest, so (ip) = 
(i/j)n ~ {i>)ii- It is customary to assume a relatively broad 
bandwidth of ~ 0.2-0.5 MHz while computing the theoretical 
signal (e.g. Zaldarriaga et al. 2004). This choice is motivated 
by the fact that it roughly matches the length scales (~ 10 Mpc) 
corresponding to the synthesized beam (4') of the upcoming 
radio interferometers, such as MWA. However, the frequency 
resolution of future experiments is much narrower, ^ 1 kHz 
(by design, this is necessary to remove radio frequency inter- 
ference [RFI]). This corresponds to scales much smaller than 
other length scales in the problem. Mellema et al. (2006) have 
shown that the power spectrum on the scales we will be inter- 
ested in (i.e. on the characteristic size <; lOMpc, of the HII re- 
gions), is insensitive to the frequency resolution once it is well 
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below the above value. Furthermore, for a statistical signal, in- 
creasing the bandwidth around a given frequency does not give 
any advantage in the signal-to-noise ratio for the detection (see 
discussion below). As shown in section[3] the scales of interest 
to us are generally greater than 20Mpc and the effect of finite 
beam width of 4' is unlikely to effect the signal significantly 
at those scales. Therefore, for simplicity, and in accordance 
with our objectives in this paper, we assume, W r (r') = 5 D (r' - r) 
and B(n,n') = 8 D (h' —n). We do note, however, that the random 
placement of the resolution pixels relative to the boundaries of 
the quasar bubbles will cause a shot noise in estimates of the 
anisotropy that could become significant, and will certainly ne- 
cessitate averaging over many bubbles to detect angular struc- 
ture. Quantifying this shot noise term is beyond the scope of the 
present paper (averaging over a large number of bubbles will be 
required, in any case, for other reasons; see the discussion in §|5] 
below). 

Under this approximation, the expression for the two-point 
correlation function is simplified to 

C{r m 6)= (AT(n hri )AT(n 2 ,r 2 )) = 

= T (n)T (r 2 )[{tl>(h,nMfi2,r2))-W 2 l (6) 

Here is the angle between the line of sight h and the vector 
ri — r 2 separating the two points, and brackets denote average 
over (say) ri (as noted above, the correlation function could be 
written equivalently as a function of the angle 8 and the fre- 
quency difference f 2 -f\)- We caution the reader that, as de- 
fined above, the angle 8 is different from the angle subtended 
on the sky. The latter is typically defined in studies involv- 
ing angular correlation functions, but we find the angle relative 
to the line-of-sight, which directly follows the shapes of the 
quasar-driven bubbles, more convenient for our purposes^] 
The correlation function of tp can be expanded as 

(^(n)V(r 2 )) = <{x H (r 1 )[l+% 1 )]}{x H (r 2 )[l + <5(r 2 )]}) 

= (x H (ri)x„(r 2 )) +2(x„(ri)<5(ri)x H (r 2 )) 

+(x H (r 1 )% 1 )x H (r 2 )<5(r 2 )) (7) 

The three and four-point functions in the expression above are, 
in general, nonzero, and cannot be easily computed, as x H is not 
a Gaussian random variable (for details, see, e.g., Furlanetto 
et al. 2004). 

In what follows, we will drop the terms that contain cross- 
correlations between the density and the neutral fraction (we 
will discuss the neglect of these terms in detail below). In this 

case, 

{W>) - (V) 2 = C,&8(rm0,z) + C< x -xl (8) 

Here £gg = (<5(i"i)<5(r 2 )) is the two-point correlation function of 
the total density contrast (which we compute using the fitting 
formula for the linear spectrum in Eisenstein & Hu 1999), and 
£,xx = (-^H(ri)x H (r 2 )), to be discussed in detail below, is the cor- 
relation function owing to the inhomogeneities of the neutral 
fractionQ To make further progress, we need to make some 
assumptions about the spatial variation of the neutral fraction. 

1 We note that finite instrumental resolution may, in practice, cause a degrada- 
tion of the precision to which the anisotropy of the correlation function can be 
extraced from the data, on scales near the resolution, since the brightness tem- 
perature fluctuations will effectively be smoothed with a filter that is defined in 
the observed (r, n) coordinates, and not in (r, 9) coordinates we use to specify 
the anisotropy. 

2 Another commonly used definition of the correlation function is £ u = 
(x H (ri).r H (r2)} - (jth) 2 - We caution the reader that the definition we adopt here 
is different. 



Before we proceed to compute the anisotropy in the above 
power spectrum due to the non-spherical shapes of the indi- 
vidual HII regions, we note that a different source of anisotropy 
will already be produced by redshift space distortions due to pe- 
culiar velocities (e.g. Bharadwaj & Ali 2004; Barkana & Loeb 
2005). This effect is analogous to the redshift-space distortion 
of the matter power spectrum (Kaiser 1987), and is described 
in the linear regime by 

tgg(r n ,8,z) ^ fairuM \l + ^l3+~0*\ + 

Hss(r 1 2,2,z)(^P+^(3 2 ^P2(d) (9) 

Here (3 ~ fl® n 6 /b, we use b = 1 throughout; P2(8) is the Legendre 
function with I = 2, and 

&5(r,0,z) = J dkk 2 P(k)j (kr) (10) 

igg{r,2,z) = J dkk 2 P(k)j 2 (kr). (11) 

Here D+(z) is the growing mode of density perturbations and 
P(k) is the matter power spectrum. We have further neglected 
the fourth (and higher) moments of the correlation function, 
which make a negligible contribution (see, e.g., Hamilton 1998). 
The anisotropy we consider below, intrinsic to the correlation 
function of the neutral fraction, will be in addition to this redshift- 
space-distortion anisotropy. 

2.2. Correlation Function of the Neutral Fraction 

While the fluctuations in the over-density <5(r) are specified 
by the cosmological initial conditions, in order to describe the 
spatial variations in the neutral fraction x H (r), we need to make 
assumptions about the process of reionization. It is generally 
believed that this transition occurred by the percolation of ion- 
ized bubbles around individual sources, or clusters of sources, 
that formed in nonlinear halos. The nature of sources of ioniza- 
tion is not clear and they could either be star-forming galaxies 
or QSOs, or both. Several works have proposed simple mod- 
els for the correlation function £ VJC of the neutral fraction in this 
picture (e.g. Knox et al. 1998; Gruzinov & Hu 1998; Santos et 
al. 2003; FZH04; Zhang et al. 2007). 

We follow the simplest model described in FZH04 (adopted 
from Knox et al. 1998), which assumes that the ionizing sources 
are randomly distributed in space. In reality, the sources are 
likely to be located at the peaks of the density field and there- 
fore clustered. In general, this will increase £ xl (by about an or- 
der of magnitude, for ionizing sources in dark matter halos that 
correspond to 2-3<r peaks; e.g. Santos et al. 2003), and can also 
modify its anisotropy (see discussion in §[5]below). In the case 
of randomly distributed identical sources with space density n, 
each surrounded by an ionized bubble with volume V, we have 
the joint probability for the ionized fraction x, = 1 —x H , 

k(nMi2» = (l-e-" v °) + e-" v °[l-e-' HV - V ° ) ] 2 , (12) 

where V = V (r\ - r 2 ) is the volume of the overlap between two 
ionized bubbles located at ri and r 2 . The first term in equa- 
tion (fT2l represents the probability that both points are ionized 
by the same source; the second term describes the case when 
they are ionized by two different sources. Defining p same 
I - exp(-nV ( ,) as the probability that the two points belong to 
the same ionized bubble (i.e., when the overlap region contains 
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at least one source), the correlation function for the neutral frac- 
tion, obtained from the above equation, can be conveniently re- 
written (Barkana & Loeb 2006) as 



x H (ri)x H (r 2 ) 



(13) 



( 1 Psame) 

We note that equation (fT2l assumes that the bubble bound- 
aries are "static" and not affected by overlap. A more realistic 
assumption would be that the total ionized volume is conserved, 
which would require expanding the joint boundaries of two 
bubbles that have merged. This would guarantee (1 -x H ) = nV, 
but modifying equation (fT2l) to take the bubble mergers into ac- 
count would require further assumptions and would be overly 
complicated. Alternatively, one could start with the "postulate" 
that bubbles do not overlap (i.e., counting a merged bubble as 
a single object); however, this condition would not allow the 
bubbles to be randomly distributed on small scales. 

The above ambiguity affects the correlation function on scales 
comparable to the bubble size, but we expect it would not sig- 
nificantly modify its predicted anisotropy. In any case, merg- 
ers between bubbles will likely dilute the anisotropy signal we 
discuss below, and our focus therefore will be at the earliest 
epochs of reionization, when the filling factor of the ionized 
bubbles is small. In this limit, nV <C 1, overlapping bubbles 
are rare, and the ambiguity is avoided: we can approximate 
Psame ~ nV„ = nfV « (1 -x n )f, where / = Vo/V. 

If the ionized bubbles were spherically symmetric with V = 
(4tt/3)R 3 then we would have the simple expression 



3r 
4Ri 



1 

— i 

12 



(14) 



for r<2R and zero otherwise, and where r =| ri— r 2 | (cf. equa- 
tion 18 in FZH04). Allowing the characteristic size R = R(z) 
and the mean neutral fraction x H = x H (z) to evolve with redshift 
introduces a new ambiguity in the redshift at which x s and / 
is to be evaluated in p same = (1 —x H )f. In practice, however, 
x H and R should evolve only on long time-scales, comparable 
to the Hubble time. Indeed, Barkana & Loeb (2006) discuss 
an anisotropy in the correlation function caused by this evolu- 
tion, and the finite speed of light. As discussed in the Intro- 
duction, our focus here will be on the anisotropies on smaller 
scales, comparable to the size of individual quasar bubbles. The 
anisotropy from the cosmic evolution is negligibly small on 
these scales (see Figure 3 in Barkana & Loeb 2006), and we 
therefore ignore it here and use the simple expressions, 



(1 Psame) 



and 



Ptamc = (\-x s )f(.r,K). 



(15) 



(16) 



Note that the anisotropy we discuss below will arise entirely 
from the dependence of /? S ame(i"i,r2) on "orientation" (i.e., the 
angle between ri - r 2 and the line of sight). 

Let us consider various limits of the expressions in equa- 
tions ( fT~5T > and ( [TBI : (i) for r tending to zero / tends to unity; 
(ii) as r tends to infinity /? sam e tends to zero and the correla- 
tion function approaches x\, as they should; (iii) as the average 
value of the neutral fraction tends to zero, the correlation func- 
tion vanishes, and (iv) and when the average value of the neutral 
fraction tends to unity, should approach as indeed it does. 




FIG. 1. — The apparent shape of the Stromgren sphere. The solid curve 
shows the locus of points, in the plane containing the ionizing source (O), 
from which light is observed to arrive on Earth (toward the bottom of the plot) 
simultaneously. The age of the Stromgren sphere is assumed to be t\\\\ = 10 7 yr, 
and the values of the other parameters are: clumping factor C = 5, neutral 
fraction x H = 0.5, source redshift z = 10 and ionizing photon luminosity S = 
10 57 sec -1 . For reference, the dashed curve shows a circle with radius R naK . 
The notation follows Yu (2005). 



2.3. Light Travel Delay Anisotropy of Individual HII Regions 

Assuming that all ionizing sources carve out spherically sym- 
metric Stromgren spheres^] in their rest frames (at least statisti- 
cally; we shall discuss this assumption in detail below), it can 
be shown that as light from different points on the Stromgren 
sphere is observed simultaneously by an observer on Earth, the 
sphere will appear distorted. The effect is, roughly, an elonga- 
tion along the line of sight; it is caused by the fact that photons 
observed at increasingly large impact parameters away from the 
line of sight to the source have traveled a longer path (Cen 
& Haiman 2000; White et al. 2003; Wyithe & Loeb 2004c; 
Yu 2005). As a result, at large impact parameters, we effec- 
tively observe the source at a younger age when it had a smaller 
Stromgren radius. 

Estimating the apparent distortion entails computing the lo- 
cus of all points lying on the Stromgren sphere from which the 
light is received at the same time. Expressions for this equal- 
arrival— time surface, as a function of the angle relative to the 
line of sight have been obtained in the case of a steady source 
(see Yu 2005 and references therein). We here adopt the results 
of Yu (2005), and include a somewhat different, brief derivation 
for completeness. 

Let us follow the notation of Yu (2005), and refer to the loca- 
tion of the quasar (or any generic ionizing source) as "(9" and 
to the point on the Stromgren sphere that intersects the line of 
sight to the quasar as "C". Let 8 refer to the angle between 
OC and OB, where B is some other location on the Stromgren 
sphere's surface (see Figure Q] for an illustration). The proper 

3 Technically, these are not Stromgren spheres, since they likely do not reach 
equilibrium within the source lifetime, owing to the long recombination time 
in the IGM (e.g. Shapiro & Giroux 1987); however, we follow here the wide- 
spread use of this misnomer. 
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(not comoving) length of OC is given by 

v(t')dt' 



(17) 



where v is the proper propagation velocity of the ionization 
front in the rest frame of the quasar, f, is the cosmic time at 
which the quasar switched on, and f, + tun is the cosmic time at 
which the quasar's photons reach the point C. Here tun stands 
for the age of the HII region at the time it is observed. We note 
that since photons take a finite time to travel between the source 
and the edge of the HII region, this age is larger than the age t q 
of the quasar, defined at the time it produced the photons that 
resulted in the HII region: tun = t q +R m!ix /c. Our aim is to find 
R(9), from the condition that photons traveling from the point B 
reach Earth simultaneously with photons from C. This requires 
that light leave the point B earlier than C by the difference in 
light travel times, 



At = -[R n 

c 



-R(6)cos(6)]. 



(18) 



(Note that the angle subtended by the Stromgren sphere is small, 
and we can safely assume that light rays reaching Earth from C 
and B are parallel.) When the ionization front reaches the point 
B, the age of quasar is smaller by this time delay, and R(8) will 
be smaller then R msix ; its value can be determined from the im- 
plicit relation 

rti+t m ,-At 

R(9) = / v(t')dt'. (19) 



with/?(0) = fl max . 

Given the velocity v(t) of the ionization front, R(9) can be de- 
termined numerically from this relation (similar to the implicit 
relation in eq. 3 of Yu 2005). A solution exists as long as v < c. 
For v = c, a solution exists for 9 < it/ 2 but not for larger val- 
ues of 9 which just means we can only see the front part of the 
Stromgren sphere. The solution R(8) tends to a constant at late 
times, and R(9) becomes independent of 9 as v/c tends to zero 
(as should generally be the case at late times). 




f?(radian) 

FIG. 2. — The comoving radius of the Stromgren sphere (1 +z)R(S) is shown 
as a function of 9, the angle between the line of sight and the vector connecting 
the position of the source with the point on the Stromgren sphere. The four 
curves show results for different ages of the Stromgren sphere: from bottom 
to top the ages are ?hii = {1, 1.5,2.25,3.75} X 10 7 years. The other parameters 
are the same as in Figure[T] 



The physical velocity of the growth of the Stromgren sphere 
v = dr/dt can be computed by (see e.g. White et al. 2003; 
Wyithe & Loeb 2004b): 



(S+Airrriuc) 



dr 
dt 



-Hr 



= c[S 



4tt 



r'n^Ca?, 



(20) 



Here r is the physical (not comoving) radius of the Stromgren 
sphere, S is the ionizing photon luminosity of the central source, 
H = H(z) is the Hubble parameter at redshift z, c is the speed 
of light, « H is the neutral hydrogen number density, cub is the 
case B recombination coefficient at the temperature 10 4 K, and 
C = / (n ml ) 2 . is the small-scale clumping factor of ionized 
gas inside the Stromgren sphere. This formula gives the correct 
limit of dr/dt — > c as r — > 0, i.e. the Stromgren sphere expands 
at the speed of light in the beginning. 

In Figure [2] we show the comoving radius (1 +z)R{9) for a 
set of fiducial parameters that could represent a bright quasar at 
z = 10, for four different ages of the HII region. As the figure 
shows, the anisotropy can be large (of order unity) for if the 
source is young (fHii=fewx 10 7 yr). 

As explained in the previous section, the effect of the 9— 
dependence of R on the correlation function is that f(r,R) be- 
comes a function of 9. Geometrically, equation ( TT4T > represents 
the volume of overlap between two spheres of radius R that are 
a distance r apart, divided by the volume of a single sphere, 
4ir/3R i . This has to be replaced by the new function f(r,9), 
which represents the normalized overlap volume between two 
objects whose shapes are described by equation (fT9l . Note that 
the shape has azimuthal symmetry, but in general, the overlap 
volume does not. We numerically compute this overlap vol- 
ume, with the line connecting the "centers" of the two objects 
(i.e. the positions 0\ and O2 of the sources) of length r, and 
oriented at angle 9 to the line of sight (divided, again, by the 
volume of a single object). The correlation function (eq. |6]l 
can now be computed using equations (|8l-(fl6]>, but replacing 
f{r,R) in equation (fl4l by the new f(r, 9). 

3. RESULTS 

The correlation function at a given redshift (eq. includ- 
ing its anisotropy caused by the light travel delay, depends on 
parameters the properties of the sources (such as their typical 
age, luminosity, and number density) and of the IGM (such as 
its density, mean neutral fraction, and clumping factor). For a 
given reionization history, the correlation function will be de- 
termined by averaging over the distribution of the ages and lu- 
minosities of the sources that co-exist in the IGM at any given 
redshift. 

We first consider, in § 13.11 simple toy models, in which all 
ionizing sources (at a given redshift) are identical, and propa- 
gate into a fully neutral IGM. These results are intended to illus- 
trate the level of anisotropy expected from particular sources, 
and the dependence of the anisotropy on various parameters. 

In § 13.21 we present more realistic estimates of the possible 
level of anisotropy. These results differ from those in the toy- 
model case, in that (i) we consider coexisting sources with a 
range of luminosities, (ii) we include the dilution of the anisotropy 
due to the presence of "fossil" HII regions, and (iii) we con- 
sider the scenario in which the HII regions produced by bright 
quasars expand into a medium that is already partially ionized 
by pre-existing galaxies (and/or lower luminosity, non-relativistic 
quasars). 

Further complications will be discussed in §|5]below. 
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FIG. 3. — The percentage chang^amathe correlation function, defined 
as A£/§ = K(O)-£(0)]/£(O) X 100, is plotted as a function of 6, the an- 
gle between the line of sight and the separation vector r/12 (note that this 
is different from the angle subtended on the sky that is commonly used 
for angular correlation functions). Each curve corresponds to a different 
value of r\2 (comoving separation between two points). The IGM is as- 
sumed to contain randomly distributed, identical ionized bubbles. The clump- 
ing factor (C = 5), neutral fraction (,v H = 0.5), redshift (z = 10) and ioniz- 
ing photon luminosity (S = 10 57 sec~') are the same for each curve. The 
four panels show the results for different values of the age of the Strom- 
gren sphere fun! clockwise from top left, the ages are {1,1.5,2.25,3.75} X 
10 7 years. The value of r\2 and the normalization of the temperature 
power spectrum, given by C(rn,9 = 0) in equation {6}, from bottom to 
top in each panel, respectively, are: ryi = {5, 10, 15,20}Mpc and {2.2 X 
KT 4 , 1.1 X 10~ 4 ,5.2 X 10- 5 , 1.2 X 10- 5 }K 2 ; m = {5, 15,20,31}Mpc and 
{2.6 x 10~ 4 ,6 x KT 5 ,5.2 x KT 5 ,4.4 x 10~ 6 }K 2 ; r n = {12,25, 35, 45}Mpc 
and {1.9 X 10~ 4 ,6.8x l(r 5 ,2.4x 10~ 5 ,1.6x 10~ 6 }K 2 ; r 12 = {35,50,58}Mpc 
and {5.6 X 10~ 5 , 1 x 10~ 5 , 1.8 x 10~ 6 }K 2 . For reference, the dashed curves 
in the top left panel show the anisotropy caused by redshift-space distortion 
in £s8 (eq. |9j. The sign of the redshift-space distortion anisotropy, as de- 
fined above, is everywhere negative, and the dashed curves shows the absolute 
value of this change. The length scale for the curves, from bottom to top, 
is {5, 15}Mpc and the normalization of 4?<5<s(0) is {7 X 10~ 6 ,9 X 10- 7 }K 2 . 
The anisotropies in £gg and in £. u are comparable, but the contribution of the 
density fluctuations £gg to temperature anisotropy (eq. [8} remains negligible 
for the the range of length-scales shown in this figure. 

3.1. Illustrative Toy Models 

In this section, we consider a population of identical sources, 
in order to isolate their contribution to the anisotropy in the 
more realistic models described in the next subsection. The 
most important source parameters are the luminosity and life- 
time; large anisotropies are expected only for bright and short- 
lived sources, such as luminous quasars. As we show below, 
the anisotropy from ionized bubbles around galaxies or fainter 
quasars is expected to be negligible. 
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FIG. 4. — All the parameters arjitdhei! notation are the same as Fig- 
ure |3] except the photon ionizing luminosity S = 5 X 10 56 sec~'. The four 
panels show the results for different values of the age of the Stromgren 
sphere tun; clockwise from top left, the ages are {1,1.5,2.25,3.75} X 
10 7 years. The value of ri2 and the normalization of the temperature power 
spectrum, given by C^r[2,8 = 0) in equation (6), from bottom to top in 
each panel, respectively, are: ryi = {5,15,20}Mpc and {2.1 X 10^,4 X 
10~ 5 ,5 x 10~ 5 }K 2 ; r 12 = {15,25,29}Mpc and {8 x 10~ 5 ,1.2 x 10~ 5 ,5.2 x 
10 _6 }K 2 ; r 12 = {25,31,38}Mpc and {5 X 10^,2 X 10~ 5 ,2.8 X 10 _6 }K 2 ; 
r 12 = {31,38,48}Mpcand {5 x 10" 5 , 2.4 x l(T 5 ,3.2x 10~ 6 }K 2 . 

3.1.1. Quasar-like sources 

Quasars can have high photon luminosities (~ 0.5-5 x 1 57 sec -1 ) 
and short ages (of a few x 10 7 years; see, e.g., Haiman & Cen 
2002, or Martini 2004 for a review). In Figure[3]we plot the cor- 
relation function for different choices of these parameters. The 
results from Figure[3]can be briefly summarized as follows. 

(a) Level of Anisotropy. The distortion from light travel de- 
lay can lead to substantial anisotropy in the correlation func- 
tion. The level of anisotropy is ~ 10% or higher, if the source 
lifetimes are ^ 2 x 10 7 years, on scales 5 Mpc <; r\2 & 25 
Mpc. The particular distortions in the shapes of the Strom- 
gren spheres translate to a characteristic ^-dependence of the 
anisotropy. 

(b) Dependence on Source Lifetime. For fixed source lumi- 
nosity, and fixed separation rn, the anisotropy generally dimin- 
ishes as the source lifetime gets longer. The sensitivity is very 
steep: e.g. for r\2 — 35 Mpc, the anisotropy for fnii = 3.75 x 
10 7 years is about an order of magnitude smaller than for fnn = 
1.5 x 10 7 years. The steep dependence of the anisotropy on the 
source lifetime is explained by the fact that the HII regions ini- 
tially expand rapidly at relativistic speeds, but then become in- 
creasingly non-relativistic as the source age increases. From 
equation d20b . we can infer that the transition occurs sooner for 
less luminous sources: the critical radius is r m ~ [5/(47r« HI c)] 1/2 
and as r ~ ct during the relativistic phase, the corresponding 
critical age of the source at the transition scales as f nr oc S 1 ' 2 . 
If sources have typical life-times that are independent of S (or 
at least scale less steeply than S 1 ' 2 ), then we can expect that 
bright sources with t <; t m will be in their relativistic phase, 
whereas fainter sources with t <^ t m will have non-relativistic 
and isotropic bubbles, with t m given by 

(21) 

(c) Dependence on Scale. The anisotropy shows a steady 
increase with increasing r\2 and then peaks at r\2 ~ 2R C , where 
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2R C = R(0)+R(n) represents the diameter of a single HII bubble. 
For larger separations, the anisotropy decreases with increasing 
rn- 

(d) Dependence on Source Luminosity. Figures[3]and[4]show 
the change in anisotropy when the luminosity of the quasars 
is changed by a factor of ten. A comparison of these figures 
shows that a decrease in luminosity is generally associated with 
a decrease in the anisotropy at a given length scale, for a given 
age of the QSO. This is to be expected as the HII region of the 
QSO remains relativistic for a longer period if the luminosity is 
increased Ceq.l2Tb. 

(e) Dependence on IGM Clumping Factor. We have repeated 
the calculations in Figures [3] and [4] with the clumping factor 
C varied in the expected range of 1 <S C ^ 20. The effect of 
an increased clumping factor is to slow down the expansion of 
the ionization fronts (see eq. [20b , and therefore to reduce the 
anisotropy. However, this effect is small unless the ionization 
front speed is limited by recombinations, which happens only 
at late times, and only if the clumping is near the high-end of 
the expected range (C ~ 20). Even in this case, we found that 
the reduction in the level of anisotropy is at most ~ 50%; the 
dependence on C is negligibly small for C <, 10 (however, the 
clumping of the IGM has the additional, and more significant 
effect of determining the ratio of fossil/active HII regions; see 
discussion of this in § !3.2l below). 

(e) Comparison to Redshift-Space Distortions. It is useful to 
compare the anisotropy caused by quasars to the redshift-space 
distortion induced anisotropy. The two sources of anisotropies 
are, in fact, coupled through the first term of equation |(8). To 
disentangle these two effects, we compute the correlation func- 
tion in the density perturbations £ss alone. The result is shown 
by the dashed curves in the upper left panel in Figure [3] This 
panel shows that the redshift-space distortion induced anisotropy 
in £ss is significant, and at a comparable or larger magnitude 
than the anisotropy in £ AA caused by the quasars. However, the 
anisotropy in the two-point correlation function is given by the 
term (1 + £gg)£xx m equation ([BJ. In the range of length scales 
shown in Figures [3] and [H £ss *C 1, so that the contribution of 
the density fluctuations £ss to the correlation function and its 
anisotropy remains negligible. We also note that for a range 
of scales, the sense of change of the correlation function, as 
a function of the angle to the line of sight, has the opposite 
sign compared with our signal from the finite light travel time. 
The redshift-space distortion induced anisotropy is well under- 
stood theoretically, and it would be interesting to estimate it 
directly from the HI signal. However, the presence of even a 
small amount of ionization-inhomogeneity induced anisotropy 
would suffice to mask this anisotropy. 

3.1.2. Star-forming galaxies 

In general, such sources are expected to be much less lumi- 
nous than bright quasars, and therefore equation ( l2Tb suggests 
that they can not cause any significant anisotropy. We can esti- 
mate the maximum photon luminosity of early galaxies by not- 
ing that 2er peaks collapsing at z = 10 correspond to halos with 
a total mass of ~ 10 8 M Q , containing ~ 10 7 M Q of baryons. If 
all baryons turn into stars with a Salpeter initial mass function 
(IMF), each proton yields about 4000 ionizing photons. This 
efficiency can be further boosted by a factor of up to ~ 20, 
if the stellar population consists entirely of massive, metal- 
free stars (see Haiman & Holder 2003 for further discussion 
and references). Assuming the ionizing photons are produced 



within a single burst lasting 10 7 years, we find S w 3 x 10 54 s _1 . 
From equation (fJTJ, we conclude that even under these ex- 
treme assumptions, the typical Stromgren spheres around these 
sources will be expanding non-relativistically. We have ex- 
plored a range of luminosities 10 53 -10 54 s _1 and ages of 5 x 
10 6 -5 x 10 7 years that could be expected for the brightest star- 
forming galaxies, and a wide range of parameters of the IGM, 
and indeed, we found that the anisotropy of the correlation 
function is dominated by the redshift-space distortion. There- 
fore if the star-forming galaxies are the exclusive agents of the 
reionization process, the light travel time delay anisotropy is 
unlikely to be detected. Conversely, this implies that any se- 
cure detection of the anisotropy would immediately reveal the 
presence of bright, quasar-blown bubbles. 

3.2. Possible Anisotropy in More Realistic Scenarios 

As stated above, a more realistic model would have to in- 
clude averaging over a range of luminosities and ages at a given 
redshift. The simple models shown in Figures[3]and[4] however, 
already allow us to draw some basic conclusions. First, since 
the anisotropy decreases steeply with source lifetime, if all of 
the ionizing sources have ages far exceeding a few x 10 7 years, 
and/or they have ionizing luminosities significantly below a few 
xl0 56 sec _1 , then the correlation function anisotropy will be 
negligible (<C 1%). This critical value for the lifetime coin- 
cides with expectations for bright quasars (see the review by 
Martini 2004), whereas the critical luminosity is about an order 
of magnitude below those of the already known SDSS quasars 
atz > 6. 

Quasars of the required brightness almost certainly appeared 
only in the late stages of reionization. These relevant later parts 
of the reionization history likely involved a mixture of both 
types of sources discussed above. The process of reionization 
was likely initiated by low-luminosity sources, such as star- 
forming galaxies (with masses £ 10 8 M Q ), or smaller black 
holes (with masses < 10 6 M Q ; Madau et al. 2004; Ricotti & 
Ostriker 2004), at high redshift (z <; 10). However, at some- 
what lower redshifts, but well before reionization is completed, 
bright quasars could have formed in higher mass dark matter 
halos (~ 10 12 M Q ) corresponding to the rarer high-(T peaks of 
the underlying density distribution. The unresolved soft X-ray 
background puts a tight limit on the contribution of quasars to 
reionization at z ~ 6 - 20, but does not rule out the above sce- 
nario, in which quasars contribute a few ionizing photons per H 
atom at z 6 (see Fig. 1 in Dijkstra, Haiman & Loeb 2004). 

The toy models above suggest that if bright quasars con- 
tribute a non-negligible fraction of the ionizing photons pro- 
duced at z > 6, then the anisotropy of their Stromgren spheres 
could be observable via the measurements of the 21cm correla- 
tion function. In the rest of this section, we present somewhat 
more elaborate models, to quantify the anisotropy in a scenario 
in which the universe is partially ionized by faint sources (star- 
forming galaxies and/or low-luminosity quasars) and the HII 
regions of bright, relativistic QSOs, with a range of luminosi- 
ties, are expanding into this partially ionized medium. 

More specifically, we make the following modifications to 
the toy models of the previous section. 

(1) Pre-ionization by non- relativistic sources. First, we as- 
sume that at the epoch of interest, non-relativistic sources (galax- 
ies or low-luminosity quasars) had already ionized 60% of the 
IGM. This is consistent with current constraints on the neutral 
fraction around z ~ 6 (e.g. Fan, Carilli & Keating 2006). The 
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size of the HII regions of these bright QSOs is likely to be larger 
than any other scale in the problem (e.g. compared to the clus- 
tering scale of HII regions of other sources) and therefore one 
can assume that at the scales comparable to the HII region of the 
bright QSO, the only source of ionization inhomogeneity is ow- 
ing to the presence of the bright QSOs. Alvarez & Abel (2007) 
and Lidz et al. (2007) present models which might bear out 
this picture (although, as these authors argue, fluctuations from 
galaxies can be important in other contexts involving smaller 
scales, such as interpreting Lyman a absorption spectra). For 
our purposes, the partial pre-ionization can be regarded as uni- 
form, and we therefore assume that its sole effect is to reduce 
the HI density in the IGM (in eq.l20b. 

(2) Additional ionization by bright quasars. The ionization 
fraction caused by bright QSOs could be a small fraction (say, 
5-10%) compared to the other sources (say, 60%, as adopted 
above). To be explicit, we set x„ = 0.05 as our fiducial choice for 
the quasar contribution, reduced by a factor of 10 from the toy- 
models. To motivate this choice, we note that each quasar with 
the smallest luminosity of interest (^ 5 x 10 56 sec _I ) would 
contribute a newly ionized volume of ~ 10 4 Mpc 3 . With the 
steepest allowed logarithmic slope for the quasar LF at z ~ 6 
(see discussion in §[5}, we expect a space density of ~ 10~ 5 Mpc 3 
for these quasars. The choice of 5-10% for the fraction of the 
volume ionized by these relativistic sources is therefore about 
the maximum allowed by the upper limits on their space den- 
sity. For reference, it is also useful to note that at z ~ 5 - 6, the 
known population of bright and detectable quasars contribute 
only about 1% of the ionizing background (e.g. Madau, Haardt 
& Rees 1999; Fan et al. 2001). Sbrinovsky & Wyithe (2007) 
used the observed quasar LF to explicitly find an upper limit of 
14.5% on the contribution of luminous (but below the detection 
threshold of SDSS) quasars at z ~ 6 to the ionizing background. 
In practical terms, this means that p sdme (eq.[T6ll and the correla- 
tion function will now be defined with respect to the ionization 
fraction caused by the bright QSOs alone (e.g. 5%). The corre- 
lation function of the HI fluctuations at the scale of the HII re- 
gions ( > lOMpc) of bright QSOs will be dominated by the ion- 
ization inhomogeneities caused by these HII regions (the other 
component will come from the density inhomogeneities which, 
as follows from Figure[3] will contribute negligibly). This also 
means that the light-travel time anisotropy at these scales will 
be mostly owing to the presence of relativistic HII regions of 
the bright QSOs, even though they contribute a small fraction 
to the ionization fraction. 

(3) A finite spread in quasar luminosities. It is obviously 
unrealistic to assume that all quasars have the same luminosity. 
To relax this assumption, we extend our analysis to include HII 
regions with a range of sizes. HII regions with volume V then 
contribute an ionized fraction nV, and one must compute the 
anisotropy by summing over the contribution from different HII 
regions. As shown in FZH04, the existence of different bubble 
sizes can easily be incorporated, by redefining quantities: 

x e = J2 x e( R ) (22) 

R 

PsameC'", 0) = same (>",/?, 9). (23) 

R 

(24) 

Here R is the radius of the HII region corresponding to a par- 
ticular type of source, determined by the luminosity and life- 
time of the objects (note that this expression is strictly valid in 



the limit we consider here, i.e. riV <C 1). For simplicity, we 
consider here only a finite spread of luminosities, in the range 
2 x 10 56 sec _1 <, S < 2 x 10 57 sec _I . Quasars brighter than this 
range are already known to be rare. There could, of course, 
be additional quasars below this luminosity range, but we as- 
sume their HII regions would be nonrelativistic, and implicitly 
lump any such quasars together with the galaxy pre-ionization 
discussed above. We keep the lifetime fixed at ?hh = 3 x 10 7 
yr (in principle, one should also allow for a spread in ages be- 
tween < fen < 3.75 x 10 7 yr, which would slightly increase 
the anisotropy). Within the 2 x 10 56 sec _I ^S^ 2 x 10 57 sec _1 
range, we consider either a flat or a steep luminosity function, 
n(L)V(L) tx L° or n(K)V(L) cx L~ 3 . The latter slope is motivated 
by current empirical limits on the slope of the quasar LF at z ~ 6 
(e.g. Fan, Carilli & Keating 2006), while the former (flat) slope 
is motivated by the possibility that the high-z quasar LF turns 
over and flattens not too far below the SDSS detection thresh- 
old. 

(4) Fossil HII regions. An important effect that will decrease 
the anisotropy we predict for short-lived quasars is the pres- 
ence of "fossil" HII regions around dead quasars. At redshift 
z = 10, for our assumed fli,h 2 , the recombination time in the 
IGM is w 10 8 (C/5) _1 years. This is longer by a factor 5 (f H n/2 x 
10 7 yr) _1 (C/5) _I than the ages of the active HII regions. The im- 
plication is that for fen = 2 x 10 7 yr and for C = 5, at any given 
time, there would be ~ 5 fossil HII regions, which have not yet 
fully recombined, for every actively ionized Stromgren sphere. 
Since these fossil HII regions are not expanding, they will not 
produce an apparent anisotropy, and they will dilute the pre- 
dicted level of anisotropy by a factor of ~ 5. Several authors 
have attempted to compute the gas clumping factor C(z) from 
first principles, using numerical simulations or semi-analytic 
arguments. The overall clumping factor at high z is dominated 
by dense gas in collapsed halos (Haiman et al. 2001) and the re- 
sulting small-scale clumping is poorly understood. Iliev, Scan- 
napieco & Shapiro (2005) find in a simulation that the contribu- 
tion from the low-density IGM to the clumping at z = 10 gives 
C = 8; the dense gas inside halos can, however, increase this 
value by a factor of a few, so that the clumping factor may 
be as high as C = 20 (Haiman et al. 2001; Iliev et al. 2005). 
On the other hand, the time taken by the fossil gas to recom- 
bine will actually vary with the density, implying that the dense 
parts of the fossil will recombine faster, while less dense parts 
will linger longer. Furthermore, the precise value of the rele- 
vant clumping factor will also be modified due to clustering of 
gas near quasars (one expects local overdensities, at least in the 
central regions of the fossil), and will also depend on the instru- 
mental specifications. The active/fossil bubble abundance ratio 
will also be modified if the formation of quasars, within the 
last recombination time, is skewed to more recent epochs. In 
our fiducial model, we assume C = 7, which, for our choice of 
?hii = 3 x 10 7 yr _I , implies that fossils outnumber active bubbles 
by a factor of ^three. We include the presence of these fossils 
in our calculations explicitly, treating them exactly as active HII 
regions, except that we assume that they are isotropic (i.e. set- 
ting p sam e to be isotropic). Note that the clumping factor, for 
the purpose of computing the speed of the ionization front (in 
eq. [20li is different, and is still taken to be C = 10. 
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FIG. 5. — The expected level of attiSptetipy)in our fiducial "realistic" model 
at redshift z = 7. In this scenario, 60% of the IGM volume is assumed to 
be already ionized by non-relativistic sources. Quasars with a range of ion- 
izing luminosities (2 X 10 56 sec~' < S < 2 X 10 57 sec~') produce HII regions 
that ionize an additional 5% of the IGM volume; isotropic fossil HII re- 
gions, left behind by dead quasars, are assumed to fill an additional 15%. 
The quasar lifetime is assumed to be fixed at fun = 3 X 10 7 yr. The clump- 
ing factor, for the purposes of the propagation of expanding ionization fronts, 
is taken to be C = 10. The scales in both panels, from bottom to top, are 
ri2 = {35,42,50}Mpc. Left panel: The curves assume a steep quasar lu- 
minosity function, corresponding to nV ex Zr 3 , and the normalizations are 
{3 X HT 6 , 1.4 X Kr 6 , 3.7 X Kr 7 }K 2 . Right panel: The curves all correspond 
to a flat luminosity function, i.e. a uniform distribution nV oc L , and the nor- 
malizations are {2.9 x lO" 6 , 1.4 x 10- 6 ,5 x 10~ 7 }K 2 . 

The results of the calculations that take into account these 
modifications are shown in Figure|5] The modifications (l)-(4) 
introduce several competing effects relative to our toy models. 
One the one hand, there is an increase in anisotropy owing to the 
more rapid expansion of the bubbles in an already pre-reionized 
medium. On the other hand, there is a decrease in the over- 
all signal coming from these bubbles, as they contribute only 
a small fraction of the total ionized fraction, and this contribu- 
tion is further diluted by the presence of fossil regions. We find 
that the effects of the finite spread of luminosities is less impor- 
tant, although the steep LF within the luminosity range consid- 
ered produces a factor of sstwo lower anisotropy (because the 
fainter, less relativistic sources will have a higher contribution). 
Overall, we find that in the range of scales we considered, the 
level of the anisotropy can still reach above ~ 10 percent. The 
detectability of this level of anisotropy will be discussed in the 
next section. 

4. NOISE AND DETECTABILITY 

The noise characteristics of interferometric experiments for 
the detection of surface brightness fluctuations have been dis- 
cussed in detail in Fourier space in different contexts by var- 
ious authors (e.g. Bowman, Morales & Hewitt 2006; Zaldar- 
riaga et al. 2004; Bharadwaj & Sethi 2001; White et al. 1999; 
McQuinn et al. 2006). The anisotropy computation lends it- 
self more readily to interpretation in real space, we therefore 
estimate here the noise characteristics for on-going and future 
interferometric experiments in real space. The surface bright- 
ness sensitivity ATb for detecting an extended source (covering 
a solid angle equal to or larger than the synthesized beam) is: 

AT B = AT A ^f- (25) 
lis 

Here flp and ils are the solid angles of the primary and the 
synthesized beams, respectively. AT&, the antenna temperature 



Here 7^ is the system temperature, Av is the channel width for 
the measurement, Af is the total integration time, and N is the 
total number of antennas. As the noise computation is valid 
for both filled antennas and dipole arrays, we refine our defini- 
tions relevant to on-going and future experiments. The primary 
beam for a filled-aperture experiment, such as GMRT, is the 
usual "field of view". For an interferometer such as MWA, the 
primary beam is the beam of each "tile" that contains sixteen 
dipoles, and therefore is the total number of such tiles. An- 
other assumption we shall make below is that the noise in each 
synthesized beam is uncorrelated. This requires uniform cover- 
age in the visibility space. Both GMRT and MWA contain long 
baselines that can give a small synthesized beam, but the noise 
in each pixel (synthesized beam) will then be correlated. There- 
fore, in both cases, the the noise computation here is applicable 
only to "antennas" within baselines £ 1 km. 

The total number of independent (i.e. uncorrelated) pixels 
for a frequency channel is A 7 ™ ~ f2p/fi,y. Assuming N u fre- 
quency channels, the total number of independent pixels in the 
data cube is N, ~ N m N v . The quantity we wish to estimate is 
the two point correlation function of the brightness temperature 
fluctuations, 

Ci 2 = <A2i(n)Ar B (r a )). (27) 

The observed temperature fluctuation ATg(ri) contains contri- 
butions from both signal and noise, but since noise is uncorre- 
lated between two pixels, for noise, £12 = 0, and therefore the 
measured £12 gives an unbiased estimator of the signal. We fur- 
ther assume that for the measured A7g is dominated by noise 
(or in other words, the signal cannot be directly imaged in the 
given integration time; see discussion in § sec:discussion), and 
in our computations below, we will equate A7g with the noise. 
The quantity of interest is the variance of the two-point function 
(eq. |27| |. This variance can be computed by using the fact that 
each pair of measurements yields an uncorrelated random vari- 
able. The variance in the average of n such pairs is the variance 
of each random variable, divided by the number of pairs (see 
e.g. Papoulis 1984), 

_2 



(28) 



n(m,0) 

Here = (A7g) is the variance of a single pair, and n(r\2,6) 
is the number of pairs for a fixed | i"i2 1 ■ Note that we explic- 
itly write n as a function of r\2 and the angle 9, to take into 
account all possible correlations in the three-dimensional cube. 
The total number of pairs for all the pixels in the data cube, 
A 7 ,, is N ttP air = N t (N t — l)/2, which allows us to write n(r\2, 9) = 
Nt, P air x f(fi2,ff), where f(ri2,9) is the fraction of pairs within 
a given range of r\2 ± Arn/2 and 9±A9/2. This gives: 

^ 2 -N vf (n 2 ,e^{ns)- (29) 

The typical value of f(r\2,&) can be roughly estimated by the 
following arguments. By definition, f(r\2,&) is the fraction of 
pairs with a given r\2 and 8. In practice, one is likely to use, 
at least in the initial stages of an experiment with small inte- 
gration times, a small number of broad bins centered around a 
set of values of r\2 and 6, to estimate the signal. For example, 
one may choose to divide the range of lengths available with a 
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particular instrumental configuration into 10 bins with roughly 
equal width Ar\2, and likewise divide the ir/2 angular range 
into 10 bins of equal A9. Using 100 such bands, a typical value 
of f(r\z,S) will then be 10~ 2 . Of course, in reality, f{r\2,9) 
will depend on r\2 and 9, but this dependence can only be com- 
puted once the instrumental configurations, and the choices of 
the bins, are specified. 

Here we list again the assumptions we have made to derive 
the above expression. The main assumption is that the noise is 
uncorrelated between pixels, or more specifically, the entire pri- 
mary beam is assumed to be filled with uncorrelated "pixels" of 
synthesized beam. This, as mentioned above, requires uniform 
sampling in the visibility plane with a small enough "grid" to 
sample the entire primary beam. Whether this can be achieved 
depends on the array configuration. As mentioned above, the 
on-going telescope GMRT, or the up-coming MWA might be 
able to achieve such uniform sampling for baselines £ 1 km (It 
should be noted that according to present MWA strategy, the 
baseline distribution is not expected to be uniform, as we im- 
plicitly assume for our estimates, but rather weighed towards 
smaller baselines (roughly a 1/r 2 distribution with a core of 
10m (Bowman et al. 2006). McQuinn et al. (2006) show that 
the noise level expected for the weighed distribution are lower 
than for the uniform distribution of baselines (Figure 6 of Mc- 
Quinn et al. 2006). Therefore our calculations here give a slight 
overestimate of the expected noise levels.) However, the statis- 
tical homogeneity and isotropy of the reionization signal will 
allow one to average over correlation function measurements 
from different primary beams (i.e. over different patches of the 
sky). As a result, even when pixels in a single primary beam 
contain correlations, it is possible to increase the number of un- 
correlated correlation function measurements, and therefore ob- 
tain the requisite number of uncorrelated pixels. On-going and 
future experiments, such as PAST (Pen, Wu, Peterson 2004) 
and LOFAR (see www .lofar.org), use a large number of 
dipole antennas (~ 10 4 ), and our noise estimate can be applied 
to these missions, as well. 

We now give here numerical estimates of the noise, using 
equation d29l ), for GMRT and MWA; similar estimates can be 
made for missions such as PAST and LOFAR. We assume the 
following parameters common to GMRT and MWA: system 
temperature T s =440K (this system temperature corresponds 
to z — 8 [Bowman et al. 2006]; we use this value throughout), 
total bandwidth, N^Av = 8 MHz (note that the expected sensi- 
tivity depends only on the total bandwidth and not on the chan- 
nel width), total observing time, Af = 10 6 sec, and f(r\2,9) = 
10~ 2 . For GMRT, we adopt a primary beam of ~ 4°, a syn- 
thesized beam of ~ 4' and assume 15 antennas, which gives 
AT B ~ 0.35 K and <5£ ~ 2 x 10" 5 K 2 . For MWA, we assume a 
primary beam of 20°, a synthesized beam of 4' and 500 anten- 
nas (or tiles), which gives A7g ~ 0.25 K and <5£ ~ 2.5 x 10" 7 K2 
(Bowman et al. 2006). 

Comparing these expected noise levels with the signal strength 
shown in Figure|5] we conclude that while GMRT can not de- 
tect the correlation function anisotropics, MWA has the capa- 
bility to do so. For instance, the expected signal at r\i ~ 35Mpc 
is larger than the expected noise levels of MWA by a factor of 
nearly eight. The detection of the anisotropy at the level of 
2-10% will require noise levels smaller by a further factor of 
nearly 10 to 50, which would require integration times between 
0.2-2 x 10 7 sec. These expectations can be considerably im- 
proved by the using expected noise levels of MWA5000 and 



SKA (Bowman et al. 2006, McQuinn et al. 2006). 

One can compare the degree of difficulty involved in detect- 
ing the anisotropy associated with finite-light-speed effects with 
the anisotropy resulting from redshift-space distortion. As seen 
in Figure [3] even though the redshift-space anisotropy can be 
large, it is generally associated with the sub-dominant com- 
ponent (density fluctuations) of the overall signal. At scales 
~ 35Mpc, density perturbation-induced signal is more than or- 
der of magnitude smaller than the signal in Figure [5] The sig- 
nal owing to density perturbations increases at smaller scales. 
However, the forthcoming interferometers LOFAR and MWA 
are better suited for detecting the signal above scales ~ 5 Mpc 
(see e.g. Figure 6 in McQuinn et al. 2006); a future instrument 
such as SKA, which is likely to have much better sensitivity at 
smaller scales, is needed for detecting the redshift-space dis- 
tortion (McQuinn et al. 2006). The anisotropy we discuss in 
this paper, however, is dominant at large scales and, in particu- 
lar, MWA is ideally suited for detecting such a signal as it has 
much better sensitivity at large scales, owing to its large pri- 
mary beam. 

The S/N estimate above asks whether the instrument can mea- 
sure the anisotropy on a particular scale r\i and in a particular 
direction 9. This estimate may be significantly more demand- 
ing than a mere detection, or a crude characterization, of the 
anisotropy. In reality, one can define much cruder measures of 
anisotropy to look for - say, dividing the angular range into a 
few "quadrants", and combining the independent power spec- 
trum measurements within each quadrant, and for several dif- 
ferent r\z. The effective number of independent (r, 9) combina- 
tions that one can combine will depend on the actual observa- 
tional strategy, and on the band power used, and we leave more 
precise estimates to future work. 

5. DISCUSSION 

In deriving the anisotropy of the correlation function due to 
light travel time delay, we assumed that the HII regions are 
spherically symmetric around a given source of UV photons. 
For many reasons, the HII regions will deviate from a spheri- 
cal shape: (a) density inhomogeneities around the source, (b) 
anisotropic emission from the source, (c) mergers of HII re- 
gions. We briefly discuss each of the these in detail below. 

Density inhomogeneities: From equation d20j >, it follows that 
during the initial ultra-relativistic phase of expansion, the ve- 
locity of the HII front is independent of density. As discussed 
above, the HII front makes the transition from the relativistic 
to the non-relativistic phase at the radius r oc n~l(x). As a re- 
sult, the HII region will expand more slowly into higher density 
regions, and develops an intrinsic anisotropy. However, as the 
density field constitutes a homogeneous and isotropic random 
process, averaging over a large number of ionizing centers - as 
will be required to measure the correlation function - will tend 
to cancel this effect. Recombinations and radiative transfer ef- 
fects will cause further intrinsic anisotropies owing to density 
inhomogeneities (e.g., Bolton & Haehnelt 2007, Maselli et al. 
2007), but these do not become important until much later, well 
into the non-relativistic phase, when the light travel time delay 
anisotropy is anyway negligible. 

Anisotropic emission: Arguments similar to the previous case 
apply equally to beamed or otherwise anisotropic emission, as 
well: any intrinsic anisotropy in the source emission will be 
uncorrelated with its orientation relative to the line-of-sight. As 
a result, averaging over a volume containing a large number of 
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ionizing centers will diminish this effect (producing no signal 
in the limiting case of an very large survey volume). 

Anisotropy owing to bubble mergers: Even if individual HII 
regions are spherical, their mergers will result in strong asym- 
metries, at least on scales comparable to the inter-distance be- 
tween the ionizing sources in a single bubble. The intrinsic 
shape after the mergers of HII regions, or the level of its anisotropy, 
is difficult to assess analytically, and would require numerical 
simulations, which is beyond the scope of our work. However, 
any anisotropy caused by mergers also constitutes a homoge- 
neous and isotropic random process, uncorrelated with the line- 
of-sight. This means that this anisotropy will be diluted as well, 
if measurements average over many clusters of ionizing centers. 

Mergers, however, will have some additional implications. 
First, the mergers will increase the characteristic bubble size, 
and therefore also the length scale where the correlations (and 
their anisotropy) peaks (e.g. FZH04). Second, the light travel 
time delay anisotropy will depend on the typical degree of syn- 
chronization between the sources in a single super-bubble. Pre- 
dicting this synchronization would require new assumptions about 
the source population, and additional modeling, but the lim- 
iting cases are easily envisioned. If all sources in a bubble 
turned on simultaneously, this would be similar to a single, 
but more luminous and anisotropic source, which will boost 
the predicted signal (additional travel time delays between the 
actual sources will again average out, if the sources are isotrop- 
ically distributed within the super-bubble). On the other hand, 
if the sources are turned on in a perfect "sequence", one after 
the other, then the effect would be similar to increasing the typ- 
ical source lifetime (by a factor that equals the average number 
of sources in a single bubble), without changing the luminosity. 
This will diminish the anisotropy signal, as discussed above. 
In particular, if the source(s) in a bubble maintain a luminosity 
similar to that of a single bright source for more than a~5x 10 7 
years, the signal will become undetectable. Whether this occurs 
would be interesting to work out in specific models for bubble 
growth. The situation should be possible to avoid when the fill- 
ing factor of ionized bubbles is small. 




log (M halo /M a ) 

FIG. 6. — The number of dark matter halos per unit redshift and solid angle, 
at two different redshifts (z = 6 and z = 10). 

In summary, the anisotropy of the correlation function does 



not depend crucially on the assumption of the sphericity of the 
HII region around a source. Detailed changes in the correlation 
function and its anisotropy, owing to merging of sources are 
hard to assess but will generally lead to an increase in correla- 
tion length. It can also suppress the signal significantly if more 
than a few bright quasar-like sources turn on, synchronized to 
within At ps t q (the lifetime of a single source), and cluster to- 
gether in a single bubble whose expansion is maintained for 

5 x 10 7 years. On the other hand, if the synchronization is ei- 
ther Af <C t q or Af 3> tq, then we expect the anisotropy to either 
increase, or remain unaffected, respectively. 

In our models, we also implicitly assumed that the edges of 
the HII regions are sharp. In reality, the edges can be blurred 
if quasar spectra are hard (Zaroubi & Silk 2005; although see 
Kramer et al. 2007 for a different conclusion) or if galaxies 
that are clustered around the quasar contribute significant ion- 
izing flux spread spatially over an extended region (Wyithe & 
Loeb 2006). In practice, however, as long as the HII region 
boundaries are sharp enough to be well defined, such blurring 
should not significantly affect our conclusions. We also note 
that at least one of the z > 6 quasars appears to be surrounded 
by a sharp HII region, with a thickness < 1 (proper) Mpc for 
a radius of ~ 5 Mpc (Mesinger & Haiman 2004). On the other 
hand, the random placement of the instrumental resolution ele- 
ments relative to the boundaries of the quasar bubbles will ef- 
fective further blur the edges of quasar bubbles, and will also 
add shot noise to estimates of the anisotropy. 

Whether the various intrinsic anisotropics can be averaged 
out depends primarily in the number of bright sources within 
the survey volume. Roughly, one expects that in order to detect 
a <5£/£ = 10% "systematic" anisotropy, in the presence of ran- 
dom intrinsic anisotropics of order unity, the survey has to con- 
tain at least > (<5£/£)~ 2 = 100 bubbles. The number of bubbles 
is highly uncertain, but we note that the lowest photon lumi- 
nosity of interest ( <; 2 x 10 56 s _1 ) corresponds to a luminosity 
that is approximately a factor of 10-50 below that of the bright 
quasars detected in the SDSS (e.g. Wyithe & Loeb 2004b). The 
slope of the quasar luminosity function at z > 6 is expected to 
be steep, but poorly constrained observationally. Richards et al. 
(2006) placed an upper limit on the slope at 4 < z < 5.4 from 
the (lack of) gravitational lensing, and found it to be flatter than 
dlogn/dlogL <; -4. This limit, applied to z ~ 6, would allow 
the existence of up to ~ 100 deg~ 2 quasars that are sufficiently 
bright to cause detectable anisotropy, or up to 3 x 10 4 sources 
in one of MWA's 300deg 2 primary beam. 

To provide another rough estimate of the possible number of 
sufficiently bright sources, in Figure [6] we show the number 
of dark matter halos per unit redshift and solid angle, at two 
different redshifts (z = 6 and z= 10), using the fitting formula 
from Jenkins et al. (2001). The abundance of the SDSS quasars 
implies that they have host halos masses of ~ 10 12,8 M Q (e.g. 
Haiman & Loeb 2001). Figure|6]shows that the slope of the halo 
mass function at these high masses is dlogN /dlogM ~ -4, 
similar (only slightly flatter) than the upper limit on the slope of 
the quasar LF. Applying this scaling, we would expect ~ 10 4 3 
sources in the 300 deg 2 area of one MWA primary beam. Be- 
cause of the steepness of the halo mass function, however, this 
estimate is very sensitive to the Mbh -Mh a i relation at high red- 
shift. For example, if the quasar luminosity scales linearly with 
BH mass L q oc Mbh, but we had the steeper scaling Mbh <x M^ Xo 
inferred for inactive galaxies at lower redshifts (e.g. Ferrara 
2002), the expected number of sources would be reduced to 
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~ 10 23 . This would still allow the detection of ~few per- 
cent anisotropy, in the presence of order unity variations in the 
shapes of individual bubbles, but only marginally. 

In our analysis above, we assumed, for simplicity, that the 
ionizing sources are randomly distributed in space. In reality, 
the sources are likely to be located at the peaks of the density 
field and therefore clustered. As mentioned above, this will 
increase the contribution of bright quasars to £„ . The cluster- 
ing of the sources, and the corresponding magnitude of the in- 
crease, could be explicitly computed in a refined version of our 
model (see, e.g., equation 19 and related discussion in FZH04). 
Here we simply note that the increase for rare sources should 
roughly trace their linear halo bias (Sheth & Tormen 1999) and 
we expect it to be about an order of magnitude for sources re- 
siding in dark matter halos that correspond to 2-3a peaks (for 
explicit calculations of the impact of source bias on the corre- 
lation function, see, e.g., Santos et al. 2003). More interesting 
in the present context is that source clustering can also modify 
the apparent anisotropy. The three-dimensional spatial corre- 
lation function of the sources themselves will appear isotropic 
(neglecting peculiar velocities). However, the probability that a 
point is ionized is given by an integral of the space density of 
ionizing sources, around the given point, over a volume whose 
shape is anisotropic (given by eq. [I9t , The probability distri- 
bution of sources within this volume will depend on both the 
distance, and the angle with respect to the line of sight, from an- 
other ionizing source. This effect may not be negligible, since 
the correlation length of the quasar distribution may approach 
the typical sizes of ionized bubbles (quasars at 3.5 < z < 5.4 
already appear to have a correlation length as large as 25 Mpc; 
Shen et al. 2007), and will have to be included in future work 
and in analyzing actual 21cm data. 

As mentioned above, in our analysis, we neglected terms 
containing the cross-correlation between density and the ion- 
ized fraction (eq.|6]l. Furlanetto et al. (2005) argued that these 
terms are likely to be small compared to the other two terms 
that we retain, £, xx £,8s and (^ xx -x^) (in eq. [8}. This is expected 
on the scale in question (~ 20 Mpc) as the leading terms we 
retain are of order 0(x H ) or 0(x 2 ) depending on the value of 
f{r,R). The density-ionization fraction cross-correlation terms 
are of order 0{Sx H ), which is expected to be smaller than the 
terms retained as 5 <C 1 and 5/x H <C 1 at scales of interest. The 
exact ratio of the cross-terms to the retained terms can never- 
theless only be reliably calculated with numerical simulations. 
In practice, within our model, (£«— x%) dominates over £, X x£,88- 
In their updated bubble-growth model, McQuinn et al. (2005) 
explicitly compared (^.u -x 2 ,) and over a range of redshifts 
and length-scales. They find that the cross-term is indeed sub- 
dominant, although only by a factor of ~ two at early stages of 
reionization and on small scales (see, e.g., their Figure 2). We 
conclude that the cross-terms will not affect our estimates by 
more than a factor of ^two, but they will have to be included in 
a more careful analysis of actual data. 

In the context of the light travel time delay anisotropy, retain- 
ing the density-ionized fraction cross-correlation can give rise 
to a new anisotropy, owing to correlation between the redshift- 
space distortion and the anisotropy due to light travel time de- 
lay. In linear theory, the redshift-space distortion can readily 
be expanded into moments of a Legendre transform, with only 
three of these moments non-vanishing (see e.g. Hamilton 1998, 
eq. (O). The light travel time delay anisotropy can similarly be 
expanded into moments of Legendre transform (this series will 



in general have a larger number of non-vanishing moments) and 
the correlation with the redshift-space distortion could be com- 
puted. 

Finally, we note that the bright quasars producing the anisotropy 
in the power spectrum should be directly detectable with a sen- 
sitive future instrument, such as the James Webb Space Tele- 
scope (e.g. Haiman & Loeb 1998). Indeed, the 21cm maps 
may help identity such quasars to begin with. With the help 
of such identification, an alternative method to search for the 
light travel time delay anisotropy would be to stack the noisy 
21cm tomographic images centered around these quasars. This 
method will require additional modeling, in order to re-scale 
the sizes of the HII regions before they are stacked. Also, if 
the emission of quasars is anisotropic, then this can introduce a 
selection effect: the optical selection will preferentially detect 
those QSOs that appear brighter along the line of sigh toward 
us. This selection effect, if unaccounted for, will mimic the ef- 
fect of the anisotropy, since the transverse directions around the 
quasar may see systematically lower fluxes. We leave a more 
detailed discussion of this stacking approach to future work. 

6. CONCLUSIONS 

The time delay caused by finite light travel time across cos- 
mological HII regions distort their apparent shapes. This effect 
may be detectable in future redshifted 21cm observations for 
bright ionizing sources during the various stages of reioniza- 
tion, and yield constraints on the luminosity and ages of the 
sources, and the neutral hydrogen density distribution in their 
surroundings (Wyithe & Loeb 2004c; Yu 2005). In principle, 
the distortion could be measured directly in tomographic im- 
ages of individual HII regions (Wyithe, Loeb & Barnes 2005). 
Direct imaging of hundreds of sources at sufficient S/N will, 
however, be unlikely to be achieved in forthcoming experiments, 
such as LOFAR, PAST, or MWA, may have to await the con- 
struction of SKA. 

In this paper, we considered the detectability of this effect 
statistically, through measuring the anisotropy in the three di- 
mensional 21cm power spectrum on a range of scales. We 
found that the anisotropy is largest when HII regions expand 
at relativistic speeds. Our results indicate that if bright quasars 
contributed significantly (i.e. around 10 % of the ionized frac- 
tion) to some stage of reionization, then the finite-light-speed 
effect could be observable in the anisotropy of the correlation 
function of the HI distribution. For quasar luminosities <; 5 x 
10 56 sec _1 and ages ^4x 10 7 years, we expect an anisotropy 
of J> 10% in the correlation function, as shown in our results in 
Figures [3] and [4] We also compare this theoretical signal with 
the noise levels expected in on-going and future radio interfer- 
ometers that seek to detect this signal. We show that on-going 
missions, such as MWA, might be able to detect this effect. A 
detection of this anisotropy would shed light on the ionizing 
yield and age of the ionizing sources, and the distribution and 
small-scale clumping of neutral intergalactic gas in their vicin- 
ity. In particular, a secure detection of this anisotropy would im- 
mediately reveal the presence of a significant number of bright 
quasars. These sources should be directly detectable with a sen- 
sitive future instrument, such as the James Webb Space Tele- 
scope (e.g. Haiman & Loeb 1998), and indeed, the 21cm maps 
may help identity such quasars to begin with. With the help of 
such identification, an alternative method to search for the light 
travel time delay anisotropy would be to stack the noisy 21cm 
tomographic images centered around these quasars. 



13 



We thank Rennan Barkana, Steve Furlanetto, Avi Loeb, and 
Stuart Wyithe for useful comments on an earlier version of 
this manuscript. We also thank Matthew McQuinn and Miguel 
Morales for useful discussions related to the sensitivity of MWA, 
and the referee for comments that improved this paper. ZH ac- 
knowledges partial support by NASA through grant NNG04GI88G, 
by the NSF through grant AST-0307291, and by the Hungarian 
Ministry of Education through a Gytirgy Bekesy Fellowship. 



REFERENCES 

Alvarez, M. A. & Abel, T. 2007, MNRAS, submitted, (asttoTph/0703740 

Barkana, R. and Loeb, A. 2001, Phys. Rep., 349, 125 

Barkana, R. and Loeb, A. 2005, ApJ, 624, L65 

Barkana, R. and Loeb, A. 2006, MNRAS, 372, L43 

Becker, R. H. et al. 2001, AJ, 122, 2850 

Bharadwaj, S. & Ali, S. S. 2004, MNRAS, 352, 142 

Bharadwaj, S. & Sethi, S. K. 2001, JApA, 22, 293 

Bolton, J. S., & Haehnelt, M. G. 2007, MNRAS, 374, 493 

Bond, R., in Schaeffer, R., Silk, J., Spiro, M., & Zinn-Justin, J. 1996, ASP 

Bowman, J. D. , Morales, M. F. & Hewitt, J. N. 2006, ApJ, 638, 20 

Cen, R., & Haiman, Z. 2000, ApJ, 542, 75 

Chen, X. and Miralda-Escude, J. 2004, ApJ, 602, 1 

Chiu, W. A., Fan, X. & Ostriker, J. P 2003, ApJ, 599, 759 

Ciardi, B. & Madau, P. 2003, ApJ, 596, 1 

Dekel, A. & Silk, J. 1986, ApJ, 303, 39 

Dijkstra, M., Haiman, Z. and Loeb, A. 2004, ApJ, 613, 646 

Djorgovski, S. G., Castro, S., Stern, D., & Mahabal, A. A. 2001, ApJL,560, L5 

Efstathiou, G. & Bond, J. R. 1999, MNRAS, 304, 75 

Eisenstein, D. J., & Hu, W. 1999, ApJ, 511, 5 

Fan, X., Narayanan, V. K., Lupton, R. H., et al. 2001, AJ, 122, 2833 

Fan, X., Narayanan, V. K., Strauss, M. A., White, R. L., Becker, R. H., 

Pentericci, L., & Rix, H. 2002, AJ, 123, 1247 
Fan, X., Strauss, M. A., Becker, R. H., et al. 2006, AJ, 132, 117 
Fan, X., Carilli, C. L., & Keating, B. 2006, ARA&A, 44, 415 
Field, G. 1959, ApJ, 129, 536 
Ferrarese, L. 2002, ApJ, 578, 90 
Freedman, W. L. et al. 2001, ApJ, 553, 47 

Furlanetto, S. R., Oh, S.-P, & Briggs, F. 2006, Physics Reports, 433, 181 

Furlanetto, S. R., Zaldarriaga, M., & Hernquist, L. 2004, ApJ, 613, 1 

Gnedin, N. & Shaver, P. 2004, ApJ, 608, 61 1 

Gruzinov, A., & Hu, W. 1998, ApJ, 508, 435 

Haiman, Z., Abel, T., & Madau, P. 2001, ApJ, 551, 599 

Haiman, Z., & Cen, R. 2002, ApJ, 578, 702 

Haiman, Z. & Holder, G. 2003, ApJ, 595, 1 

Haiman, Z., Rees, M. & Loeb, A. 1997, ApJ, 476, 458 

Hamilton, A. J. S. 1998, The Evolving Universe. Selected Topics on Large- 
Scale Structure and on the Properties of Galaxies, Astrophysics and space 
science library (ASSL) Series vol no: 231, p. 185 

Iliev, I.T., Scannapieco, E. & Shapiro, PR. 2005, ApJ, 624, 491 

Jenkins, A. et al. 2001, MNRAS, 321, 372 

Kaiser, N. 1984, ApJL, 284, 9 

Kaplinghat, M., Chu, M., Haiman, Z., Holder, G. P., Knox, L. & Skordis, C. 

2003, ApJ, 583, 24 
Knox, L., Scoccimarro, R., & Dodelson, S. 1998, Phys. Rev. Lett., 81, 2004 
Kogut, A. et al, 2003, ApJS, 148, 161 
Kramer, R., et al. 2007, in preparation 

Lidz, A., McQuinn, M., Zaldarriaga, M., Hernquist, L., Dutta, S. 2007, ApJ, 

submitted, astro-ph/0703667 
Madau, P., Haardt, K, & Kees, M. J. 1999, ApJ, 514, 648 
Madau, P., Meiksin, A., & Rees, M. J. 1997, ApJ, 475, 429 
Madau, P., Rees, M.J., Volonteri, M., Haardt, E, Oh, S.P., ApJ, 604, 484 
Martini, P. 2004, in Carnegie Observatories Astrophysics Series, Vol. 1: 

Coevolution of Black Holes and Galaxies, ed. L. C. Ho (Cambridge: 

Cambridge Univ. Press), p. 169 
Maselli, A., Gallerani, S., Ferrara, A., & Choudhury, T. R. 2007, MNRAS, in 

press, astro-ph/0608209 
McQuinn, M., Furlanetto, S. R., Hernquist, L., Zahn, O., & Zaldarriaga, M. 

2005, ApJ, 630, 643 

McQuinn, M., Zahn, O., Zaldarriaga, M. Hernquist, L., & Furlanetto, S. R. 

2006, ApJ, 653, 815 

Mellema, G., Iliev, I. T., Pen, Ue-Li, Shapiro, P. R. 2006, MNRAS, 372, 679 
Mesinger, A. and Haiman, Z. 2004, ApJL, 611, 69 
Mesinger, A. and Haiman, Z. 2006, ApJ, in press, astro-ph/0610258 
Miralda-Escude, J., & Ostriker, J. P. 1990, ApJ, 35TJ7T 

Morales, M. E, Bowman, J. D., Cappallo, R., Hewitt, J. N, Lonsdale, & Colin 

J. 2006, New Astronomy, 50, 173 
Oh, S. P. 2001, ApJ, 553, 25 

Papoulis, A. 1984, Probability, Random Variables, and Stochastic Processes, 
(McGraw-Hill) 



Page, L. et al. 2006, astro-ph/0603450 

Pen, U.-L., Wu, X.-P., Peterson, J. 2004, astro-ph/0404083 

Perlmutter, S. et al. 1999, ApJ, 517, 565 

Richards, G. T., et al. 2006, AJ, 132, 967 

Ricotti, M. & Ostriker, J. P. 2004, MNRAS, 352, 547 

Riess, A. G. et al. 1998, AJ, 116, 1009 

Santos, M. G., Cooray, A., Haiman, Z., Knox, L., & Ma, C.-P. 2003, ApJ, 598, 
756 

Sbrinovsky, J. A., & Wyithe, J. S. B. 2007, MNRAS, 374, 627 

Sethi, S. K, 2005, MNRAS, 363, 818 

Shapiro, P. & Giroux, M. L. 1987, ApJ, 321, LI 07 

Shaver, P. A., Windhorst, R. A., Madau, P. & de Bruyn, A. G. 1999, A&A, 345, 
380 

Shen, Y., Strauss, M. A., Oguri, M. et al. 2007, AJ, 133, 2222 
Sheth, R. K, & Tormen, G. 1999, MNRAS, 308, 119 
Spergel, D. N. et al. 2006, ApJ, submitted, astro-ph/0603449 
Tozzi, P., Madau, P., Meiksin, A. & Rees, M. J. 2000, ApJ, 528, 597 
Tytler, D., O'Meara, J. M., Suzuki, N., & Lubin, D. 2000, Physics Reports, 
333, 409 

Venkatesan, A., Giroux, M. L. & Shull, J. M. 2001, ApJ, 563, 1 
White, M., Carlstrom, J. E., Dragovan, M., Holzapfel, W. L. 1999, ApJ, 514, 
12 

White, R. L., Becker, R. H., Fan, X. & Strauss, M. A. 2003, ApJ, 126, 1 

Wyithe, J. S. B. & Loeb, A. 2004a, Nature, 432, 194 

Wyithe, J. S. B. & Loeb, A. 2004b, Nature, 427, 815 

Wyithe, J. S. B. & Loeb, A. 2004c, ApJ, 610, 1 17 

Wyithe, J. S. B. & Loeb, A. 2006, ApJ, submitted, astroph/0607246 

Wyithe, J. S. B., Loeb, A. & Barnes, D. G. 2005, ApJ, 634, 715 

Yu, Q. 2005, ApJ, 623, 683 

Zaldarriaga, M., Furlanetto, S. R. & Hernquist, L. 2004, ApJ, 608, 622 

Zaroubi, S., & Silk, J. 2005, MNRAS, 360, 64 

Zhang, J., Hui, L., & Haiman, Z. 2007, MNRAS, 375, 324 



